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We consider one model of a black hole radiation, in which the equation of motion of a matter field is 
modified to cut off high frequency modes. The spectrum in the model has already been analytically 
derived in low frequency range, which has resulted in the Planckian distributin of the Hawking 
temperature. On the other hand, it has been numerically shown that its spectrum deviates from the 
thermal one in high frequency range. In this paper, we analytically derive the form of the deviation 
in the high frequency range. Our result can qualitatively explain the nature of the numerically 
calculated spectrum. The origin of the deviation is clarified by a simple discussion. 
PACS number(s): 04.70.Dy, 04.62.+V, 04.70.-s 



I. INTRODUCTION 

There are two important theories in modern physics: 
quantum mechanics and the theory of general relativity. 
The former has been applied to field theory in flat space- 
time and marked triumphs in understanding the wide 
range of microscopic phenomena. On the other hand, 
the latter has given us insights into phenomena having 
close relation to gravity. The theory of general relativ- 
ity, which is a classical theory, has drastically changed 
our understandings of space and time. Yet, we have 
no idea to treat the gravity in Planck scale. It is un- 
clear whether the unification of quantum mechanics and 
the theory of general relativity will be accomplished in 
the context of quantum physics. Today, however, we do 
not have any idea superior to the quantum mechanics in 
treating Planck scale physics. It is, therefore, meaningful 
to study how quantum and/or semi-classical effects will 
be expected in the presence of the strong gravity. 

The most remarkable discovery including a semi- 
classical gravitational effect is the Hawking radiation [Q , 
which is concluded by treating matter fields on spacetime 
as quantum while a black hole metric as classical. Ac- 
cording to this theory, a black hole radiates particle flux 
of a thermal spectrum, whose temperature is k/2t: where 
K is the surface gravity. 

The derivetion of the original Hawking radiation re- 
lies on the extremely high, over Planck scale, frequency 
modes. These modes arise from the extremely high grav- 
itational red shift which the radiation undergoes during 
propagating from the event horizon to the asymptotically 
flat region. It is the investigation of the Hawking radia- 
tion that sheds light on Planck scale physics. 

In the study of this physics, one model of a black hole 
radiation (|] @| |^ has been proposed, which makes an as- 
sumption on the effects of unknown Planck scale physics 



as follows: the equation of motion (EOM) of matter field 
is modified to have a cut off of the high frequency modes 
with respect to free fall observer. ^ There are two im- 
portant properties of this model. First one is that the 
modified EOM of matter field violates Lorcntz invariance. 
Secondly, the origin of the black hole radiation we con- 
sider in this paper, which is called the mode conversion, 
^ is different from that of the original Hawking radiation. 
We call the treatment of the black hole radiation in this 
paper as the mode conversion (MC) model hereafter. 

The spectrum of the fiux in two dimensional MC model 
has already been numerically calculated . It has been 
shown that, although the spectrum almost agrees with 
the Planckian distributionits, it deviates from the ther- 
mal one in the high frequency region k <ti lu, where uj is 
the energy (frequency) of a massless scalar field. On the 
other hand, an analytical calculation Q has resulted in 
the thermal spectrum, where the analysis has been car- 
ried out perturbatively up to lowest orders of k and lu in 
the range uj < k. This is appropriate because the deriva- 
tion of the spectrum in Q does not take the effects of the 
high frequency range into account. 

The purpose of this paper is to extend the analytical 
calculation in Q to the high frequency range, k < ui. We 
obtain the form of the distribution function, N{uj), as 



N{uj) 



1 - au?l2 



exp[(27rw/K)(l - awV2)] - 1 



which shows the deviation from the thermal spectrum. 
Here, a is the square of the cut off scale characteriz- 
ing the new physics in Planck scale. This result denotes 
the same tendency as what numerical calculations show. 
Furthermore, we will give a simple explanation of the 
occurrence of the deviation. 
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* This model was motivated by the dumb hole which is 
a hydrodynamical analogue of a black hole radiation. 

t The phenomenon called mode conversion has already been 
known in plasma physics. @ 
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In section ||, we explain the MC model |^ . Section [II 
is devoted to the calculation of the spectrum. Lastly, we 



summarize the results in section [V 



Throughout this paper, we use Planck units, h — c ^ 
G — 1, and restrict our discussion in two dimension for 
simplicity. 



II. MODE CONVERSION (MC) MODEL 

We consider the static spacetime with the metric of the 

form ii a a 



{dx — v{x) dt ] 



(1) 



It is obvious that the timelike Killing vector field of this 
spacetime is dt- This spacetime has an event horizon at 
Xh satisfying v{xh) = —1- The surface gravity of this 
spacetime, k, is calculated to be k = {dv / dx){xh) which 
is 1 /AM for the case of the Schwarzschild black hole of 
mass M . Note that the world line of dx — v{x) dt = is 
normal to the surface, t = const, and that it is a geodesic 
of a free fall observer. The velocity of the free fall ob- 
server with respect to the rest observer is v{x), where we 
require v(x) < and monotone increasing dxv{x) > 0, 
that is, v{x oo) — > vq where — 1 < I'd < by defini- 
tion. The proper time of this free fall observer is t of the 
coordinate (|l]). 

The action of a massless scalar field ^ @ @ is 



d^x^J~^ [(u($))* u($) - (s($))* s($) 



(2) 



where g = det , u is a unit tangent to the free fall world 
line and s is a unit space-like vector which is outward- 
pointing and orthogonal to u. We choose the deriva- 
tives of <& along u and s as u(<I>) — and s($) = 
F{s^dfj.)^, where F{s^d^) is a function of the differential 
operators, 9^, and we require F(— s'^c?^) = — i^(s^9^). 
This modification, i^(s''9^), expresses the MC model's 
assumption about the effects of unknown Planck scale 
physics. This leads the EOM of $ in the coordinate of 
(111) to be 



(3) 



Hereafter we specify the form of F{dx) as F'^{dx) — 3^ + 
ad^, where a = l/k^ and kc is the cut off scale which 
characterizes the new physics in Planck scale. We set the 
order of a is of unity, 0(a) ~ 1. 

We require that, at least far from the event horizon, 
v{x) is nearly constant enough to require the validity of 
WKB approximation. This means that dxv{x)/k(x) ^ 1 
and dxk{x)/k^{x) <C 1 (fc ~ const.) far from the event 
horizon. With setting $ ~ exp {—iojt + ikx), the EOM 
(0) gives the dispersion relation 



uj' = —vk + uj 

uj' = ±VA;2 - ak^ = F{k) 



(4) 



where uj' and lo are the frequencies with respect to the 
free fall observer and to the rest observer respectively, 
since u(<i>) — —i{uj — vk)^ — —iuj'<^ and — — 
The Killing frequency, tj, is conserved during time evolu- 
tion but the free fall frequency, w', is not. Equations (^) 
are the dispersion relation of this MC model, which have 
four mode solutions for fixed uj. We call the wave num- 
bers of these modes as fc_, /c_s, fc+s and fc+ in increasing 
order, as shown in the figure |l|. It is not the signature of 
u) but of Ld' according to which we can judge whether the 
solution is the positive frequency mode or the negative 
frequency one 




CO ' = F(k) CO ' = -V k + CO 

FIG. 1. Dispersion relation of our MC model. 



The origin of the black hole radiation in this model can 
be understood by analyzing a wave packet propagation 
B Hi • The group velocity of a wave packet with re- 
spect to the free fall observer and the rest observer are 
expressed as Vg = duj'{k)/dk and Vg = Vg + v{x) respec- 
tively. In the case of cj > 0, we can find easily that k-s,± 
modes are of ingoing, while fc+s one is of outgoing. In 
order to understand the physics of mode conversion, the 
hydrodynamical analogue of black hole radiation is help- 
ful |^ 1^. In this model, a fluid flows toward a center, on 
which we consider a wave packet propagation. Here, v{x) 
in (0) represents the velocity of the infalling fluid flow. 
As shown in the figure ^, k^s mode travels outward away 
from the event horizon, and against the infalling fiow 
of back ground fluid. In tracing the mode backward in 
time, it approaches toward the event horizon. In the fig- 
ure ^ we show how fc+s mode moves in the diagram of 
the dispersion relation. Near the event horizon, its group 
velocity, Vg'(fc+s), eventually coincides with the fluid ve- 
locity, v{x). Here the mode conversion takes place. In the 
view of our tracing backward in time, the infalling wave 
packet of fc+s never reaches the event horizon and is con- 
verted to the other wave packet of k-s,±- We should note 
that the negative free fall frequency mode of k_ comes 
to arise in the process of the mode conversion, while the 
other modes of k-s^-^- are positive ones. 
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FIG. 2. Schematical description of the mode conversion on 
an asymptotically flat spacetime. I* are the future and past 
null infinity, io is the spatial infinity. is the event horizon. 
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FIG. 3. Graphical description of the mode conversion in 
the graph of dispersion relation. In view of tracing the wave 
packet of fc+s backward in time, it is converted to the wave 
packet of fc_s,±. 

The situation we consider in the figure ^ where no 
transmitted wave acrosses the event horizon, corresponds 
just to a boundary condition for the calculation of the 
original Hawking radiation. In the asymptotically flat 
region, a; — > oo, we can decompose the wave function as 

^(x^co)^ Q exp(— zcj + /c/x) , (5) 

l=±s,± 

where each coefficient, c±s^±, can be uniquely determined 
by this boundary condition. Then we can obtain the 
form of the black hole radiation spectrum observed at the 
asymptotically flat region, N{lS)^ by Bogoliubov transfor- 
mation js). It becomes 

A^H = (0.1 4a/ |0,) 

c.'(fc+,)Fg(fc+.)4,(u;) ' ^ ^ 

where |0i) is the vacuum state of the free fall observer at 
the initial time (i^— oo), at and a/ are the creation and 



annihilation operators of free fall observer at the flnal 
time (t — > oo). In deriving (|^), we have assumed the 
Killing frequency spectrum of the wave packet, Lo(k), is 
sharply peaked around the four wave numbers A:±s,±- 

III. BH RADIATION SPECRTUM IN MC MODEL 
A. Analytical calculation of the spectrum 

1. The plan to calculate N{u}) 

We consider only the mode function corresponding to 
the wave number at the peak of u}{k). We set its form as 

x) = exp(-ia;t) (j){x) , (7) 

where uj is the peak value of the Killing frequency spec- 
trum of the wave packet. Then EOM becomes 

2v{iuj-v')(j}' + {uj^ +iujv')(j}^Q. (8) 

As mentioned at the end of previous section, we adopt the 
boundary condition of total reflection, which is expressed 
as 

B.C.: The solution of ^ dumps with decreasing x 
inside the event horizon. 

We take the origin of spatial coordinate x at the location 
of the event horizon, that is, v{xh = 0) = —1. 

To calculate N{u)), we follow the same plan as in [Q. 
This plan consists of three steps: First step is to solve 
(^) around the event horizon. We use the Laplace trans- 
formation in this step, where the contour of the Laplace 
integral should be chosen as it satisfies the B.C. Secondly, 
we look for the solution of (||) far from the event hori- 
zon by means of the WKB approximation. We match the 
WKB solutions with those of the first step. The existence 
of the overlap region (matching zone) should be checked, 
where the solutions in first and second steps are valid. 
Third step is to calculate N{uj) at the spatial infinity 
with using (^. 

2. Formulas for each step and small parameters 
We expand v{x) around the event horizon as 

v{x) ~ — 1 + , (9) 

1 — v{x)'^ ~ 2kx , 

where we require \kx\ <^ 1 and n is the surface grav- 
ity which is assumed to be small itself. We use Laplace 
transformation of 4>{x) , 

0(x) = / dse'''4>{s). (10) 
Jc 
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With approximating every coefficient in EOM (^) up to 
0{k) , we obtain 



(11) 



where /(s) = + iujs, B{s) = as^ — Kf3f'{s) and /3 ~ 
1 + iuj/K. This gives us 



0(x) 



ds (s + icj) 



' + liOS) X 



exp 



a /I 
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I — s — w s — ^-Q^ 



ds exp {h{s)) 



(12) 



We will find later that the factor (s + iw)*""^/^", which is 
ignored in [Q, will make important contributions to the 
deviation of the MC spectrum from the thermal one. 
For WKB approximation, we set 



(f>{x) — exp 



i / k{x')dx' 



(13) 



where xq is a constant. We expand k{x) = k(){x) + 
{1/X)ki{x) + ■■■ in introducing the scaling a; — > Ax 
(9a; -> (1/A)9^) 1], then we obtain from (|) that 



0(1) : a/co - (1 - «^)fco ^ 2tJufco + uj^ 
0(A-i):fci = ^^ln [V^(2afc3 







(14) 



-{1 - v^)ko ~ Lov)] , (15) 

where four solutions of k correspond to k±_±s. 

To calculate N{uj), we should evaluate w'Vg. This is 
given by (0) to be 



j'Vg = -2akl + (1 - -Uo)fco + t^-yo ■ 



(16) 



of integrand in ( |T^ ) in using the asymptotic expansion of 
X by the method of steepest descent contour (SDC). We 
can obtain the matching zone from these inequalities. It 
becomes 



a\i/3 1 
< a; < - , 



(18) 



where we use the forms of sq a-nd fco^i which are calculated 
later in appendix^ and at (p^). Furthermore, we require 
LUX <^ 1 for later convinience. 

With introducing A and ^ hy kx — (l/ny/a)^ and 
Lu — k{\/ K^/a)^ where n^Ja <C 1, the parameter region 
we discussed above can be expressed as 



3A>2(/i-l),--<A<0,0</i,A + ^<0, (19) 
o 

where first and second inequalities obtained by ( |l7| ) with 
(0) and (|l^) respectively. 



Parameter 
Region 




Matching Zone 
FIG. 4. Parameter region we consider. 

3. Calculation of N{uj) 



As stated in section |, N{uj) deviates from the ther- 
mal spectrum in the high frequency region k < [|| . So 
we may not be able to derive the deviation with expand- 
ing fco around = as done in Q|. Thus we focus our 
attention on the region, k < w, in order to derive the 
deviation. It is expected that the deviation becomes vis- 
ible at the value of to which does not exceed k so much. 
Thus, when we analyze how the MC spectrum deviates 
from the thermal one, the fourth term in the left-hand- 
side of (|l^) can be neglected in obtaining WKB solutions. 
This means to require 



a/cp > (1 - v^)kQ > uj^ , Lovko > uj^ 



(17) 



The other restrictions on the region of n and co come 
from the validity of the solution of WKB and that of the 
Laplace transformation. These restrictions are expressed 
by |fci/fco| < 1 from WKB and by kx < 1 and |sox| > 1 
from Laplace transformation, where sq is a saddle point 



We can proceed to calculate N{ui) with above prepara- 
tions. As we carry our the Laplace integral of (|lj) in the 
appendix we obtain the solution of in the matching 
zone as 



(x) ^ 01 (x) + 4>2ix) + (paix) , 



(20) 



where, as demonstrated below, 0i,2,3 should be matched 
with WKB solutions, (l)±,+s, respectively. These are ex- 
pressed as 



/^gi7r(n-hl/2-e3/4)-illQa)"'/12K ^ 



S7ra)/2K(l-Qa)^/2) [ ^ 

a 



x-3/4— /2- exp 



2 /2k o/, uj 
sz-^-x3/ +.-X 



(21) 



4 



= e-5«/3 2 sinh 



^(1 



r 



1 - 



auj 



',U) I K 



(22) 



where (e, n) = (—1,0), (1,1) for and (/)2 respectively 
and ^ = zQ!a;'^/4K. The factor 1 — aw^/2 shows the differ- 
ence form the analysis in |^ , which arises from the factor 
(s + iu;)*"'^'V2K in the integrand of (|l2|). 

Next we shift our step to obtaining WKB solutions. 
We can obtain fco from ( p^ ) with the conditions of ( p^ ) 
only. Then fco becmoes 



w/(l + w) 



±V(r^^^^^)7a + ^;w/(l - w^) for fc± 



where these forms of fco are valid in the region 3A > 
2{fi — 1) which is given by (^7|). Thus, for the purpose 
of matching with (|2(]|), w e obtain WKB solutions in the 
parameter region of (|19D as 



4>±{x) ~ A± (±2k) 



-3/4 ^-3/4-ii^/2K ^ 



exp 



±- 



.3/2 



(24) 
(25) 



where A±,+s are the coefficients obtained by in (|13|). 

By matching WKB solutions, ( p^ and (|5|), with that 
of Laplace transformation, (|20|), we obtain the wave func- 
tion as 



/^g«V4-ll?/3 X 
-7ruj/2K(l-cn^^/2) 



-l/4-iLj/2K-|-^ 



(-2k)3^ 



0_(a;) 



_^g-i57r/2+7rt^/«(l-Qa;V2) 0^ (j;) 



+ e"^«/3 2 sinh 

r 



^(1 



—i— 1 

K V 2 



1 



A, 



(26) 



Note that we can use this result beyond the matching 
zone, i.e. (ajn^^l'^ <^ x which is given by the validity of 
WKB approximaton, |fci/A:o| ^ 1- 

To obtain N{u!), we need to evaluate CLi'Vg and 0(x) 
at the asymptotically flat region. We obtain with taking 
the limit u(x — > oo) — > that 

(/>±(X OO) ~ A± (±l)-l/2(i _ w2)-3/4 X 



X exp 

b+si^ — > cx)) ~ exp 



±i 



1 - . 

a 1 ~ ^0 



. 2wow 



(27) 
(28) 



and 



Lj'(fc+s)Vg(A;+s) ^ -t^- 



(1-t;2)3/2 



(29) 
(30) 



Then we get c_ and by (26), (27) and (p8[), and finally 
obtain N{uj) by (§) to be 



7V(w) 



I - 



exp[(27rw/K) (1 - auj^/2)] - 1 



(31) 



This result shows that MC spectrum is enhanced in com- 
parison with Hawking spectrum. 



B. Comparison with numerical results 

We compare the analytically derived form of the spec- 
trum, (|3l]), with numerical calculations. The numerical 
calculations are carried out with two types of t;(a;), which 
are shown successively in the following two subsections. 

We have used Mathematica to solve (p|) and fit its 
solution with the form, = X]i=± ±s 6xp(iA:;), in order 
to get c_ and c+s. 

Note that we have set k/2ti — 0.0008 and a — 1 
througout these numerical computations. 



1. First type of v{x) - CJ type 

First example of u(x), which we call Corley-Jacobson 
(CJ) type takes the form 



v[x) — sgn(a;) -•\/tanh [(2^2;)^] 



1 



(32) 



To let fc_ mode be a negative free fall frequency one at 
the asymptotically flat region, the value of uj is restricted 
by upper value tOmax ~ 0.16. Further, in our parameter 
region (|l9|), the range of to is bounded as 0.005 < lo < 
0.04 (0 < M < 2/5). 

We have calculated the relative deviation, 



R{u) 



N{uj) 
Nh[uj) 



1, 



(33) 



where Nh is the Hawking spectrum. Figure |^ is the plot 
of this relative deviation. R{uj) grows qualitatively more 
and more with increasing uj. 
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0.01 0.02 0.03 

FIG. 5. Graph of relative deviation, R{uj), with CJ type of 
v{x). R{lo) grows qualitatively more and more with increasing 



This shows that the MC spectrum is enhanced in com- 
parison with the Hawking spectrum, as is indicated by 
the analytical result (|l|) . However figure ||, which is the 
plots of the logarithm of relative deviations, exhibit that 
the order of R{lo) in CJ type differs largely from our re- 
sult (|3^). We expect that this difference can be reduced 
by changing the form of v{x). The clue of this change is 
in the EOM (^. The most effective term for determining 
the analytic! ty of the solution of (H) is the second term, 
(1 — v'^)cj)" . [J Thus, in order to compare (31) with a 
numerical calculation, the expansion 1 — ~ 2kx near 
horizon should be better than the expansion v ~ —1 + kx. 
By the way, v(x) of CJ type can be expanded near hori- 
zon as: 



v{x) ~ -1 + Ka; + O ((kx)^) , 
l-v{xf ~ 2KX + 0{{Kxf) . 



(34) 



While CJ type is good form for the expansion v ~ —1+kx 
near horizon, however it is not good one for the expansion 
l — ^ 2kx. We should modify v{x) to let the expansion 
1 — ~ 2kx near horizon be better than CJ type. 



2. Second type of v{x) - DS type 

The modified type of v{x), which is called as Double- 
Squareroot-Tanh (DS) type hereafter, is 



v{x) = — 



1 -sgn(a;)- 



\ 







tanh 









(35) 



The range of u is 0.005 < uj < 0.04 (0 < /i < 2/51 
because of the same reason as that of C J type. Figure g 
shows the shapes of v{x) of CJ and DS types. They vary 
from —1 at a: = to X = —1/2 at x — *■ oo. 




FIG. 6. Shapes of v{x) in the form of CJ and DS types. 
They vary from —1 at x = to a; = —1/2 at a; oo. 

The v{x) of DS type takes the form near horizon: 

V{x) ~ ~\ + KX + O {{Kxf) , (36) 
l-v{xf ~ 2KX + 0{{Kxf) . 

The expansion v ~ ~1 + kx near horizon for CJ type 
is better than that for DS type, however the expansion 
1 — ~ 2kx for DS type is better than that for CJ type. 

Figure H is the plot of the relative deviation for DS 
type. R{uj) grows qualitatively more and more with in- 
creasing uj. This figure and the figure |^ implies that the 
details of the MC spectrum depend strongly on the form 
of v{x). The total behavior of the figure 0, however, 
shows that, althogh some data take negative value, the 
qualitative nature of the MC spectrum is enhancement 
in comparison with the Hawking spectrum, which is in- 
decated by the analytical result (pl|). Then the figure 
^ involves the plot of the absolute value of R{uj) in DS 
type. This shows that the order of R{uj) can be explained 
by our analytical result (^TJ) fairly good, since the higher 
order effects by 1 — in DS type are weakened to agree 
with the assumption made in our analytical treatment of 
the MC spectrum. 

From above it is valid to mention that our resultant 
spectrum (^l|) can explain the qualitative nature of the 
MC spectrum in the range, k < uj < K{K^/a)~^^^. 



** This is also implied in section IV 
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FIG. 7. Graph of relative deviation, R{lij), with DS type of 
v{x). R{uj) grows qualitatively more and more with increasing 

LO. 




0.01 0.02 0.03 0.04 

CO 

FIG. 8. Graphs of logarithm of absolute value of the rel- 
ative deviations. The solid line is for our analytical result. 
These show that our resultant spectrum explains the qualita- 
tive nature of the MC spectrum. 



IV. SUMMARY AND DISCUSSION 

We have derived the correction to the black hole radi- 
ation spectrum in the MC model in the frequecy range, 
< ^ < 2/5, w = nil/n^/ay. The MC model cut the 
high frequency mode off, so one may guess naively that 
the spectrum would be suppressed. Our result, however, 
shows that the MC spectrum is enhanced in the frequency 
range, < /i < 2/5. This is consistent with the numerical 
calculations as seen above. 

The deviation factor 1 — Q!a;^/2 in (|3l] ) is understood as 
a correction of red shift which outgoing modes undergo 
during propagating from the event horizon to the asymp- 
totically flat region. This is reflected on the form of fc+s 



mode, /2) This form can be easily 



explained. Neglecting the fourth term in (||) which is of 
second order of small quantities k and oj, we obtain 

aip'" + 2Kxip' - 2(1 - Kx){iuj - k)<^ = , (37) 

where we set (j)'{x) = ip{x) and v{x) ~ — 1 + kx. Under 
the assumption a < 1, we expand ip = ipQ + aipi + ■ ■ ■ 
which leads 

2Kx<y9Q — 2(1 — Kx){iuj — K)ipo = (38) 
2KXip'i — 2(1 — Kx){iuj — k)lpi = — ((5q" . (39) 



The lowest order solution is given by 

(^o(a;) =a;*"/«-ie('=-*")^. (40) 
We are interested in the analyticity of the following form: 



(41) 



since the temperature of the original Hawking radiation is 
determined by the form of an outgoing mode which takes 
the same form as above. Then we consider the correction 
to this analyticity. With setting B = iuj/K — l^wc obtain 
ifi as 



ipi{x) 



- 3B^ + 2B 3(^3 - S2) 



3^0;'^ 



3kB^ 



..2 n3 



B-'lnx 



2x2 
(po{x) . 



(42) 



We can recognize the fourth term in the square brackets 
as the secular one. Although this term is not dominant, 
it can change the analyticity which determines the tem- 
perature of the radiation. We use the renormalization 
group method in order to take all of the contributions 
of the secular term of ipi into account [^. In the 
parameter region lu/k,> 1, we obtain 



if . 



X 



i(c^/K)(l-QwV2) 



X 



(43) 



which shows the same deviation as in (31) 



Further, our result agrees with that of ||11[. In the 
paper |l^, the order of the MC spectrum's deviation has 
been estimated, which is consistent with our spectrum, 
N{oj). 
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APPENDIX A: CALCULATION OF LAPLACE 
INTEGRAL 



We approximate ( |12| ) in this appendix by an asymp- 
totic expansion of x with the method of steepest descent 
contour (SDC). The saddle points, so, are the solutions 
of dh/ds = : 



asg + 2k,xsq + 2 {iiUKx — k(3) sq ~ iujnP — . 
The (j){x) approximated by SDC is expressed as 

q>{x) ~ i\x\h"iso)r'^^ X 



(Al) 



exp 



xh{so) 



"7 + n 



(A2) 



where n is the integer determining the direction of the 
SDC and 7 = 0, 1 for x < and a; > respectively. 

With approximating (Al) and (A2) in the paramater 
region (|l^), we obtain the saddle point satisfying B.C. in 
a; < as 



So — \/2K\x\/a — i(3/2\x\ , 
and for a; > 0, the saddle points are 



So 



±i^/2Kx/a~ l3/2x 



(A3) 



(A4) 



The integral contour is determined as shown in the figure 
^, where the contour should asymtpote to one of the three 
regions, n/6 < arg(s) < 7r/2, 57r/6 < arg(s) < 77r/6 and 
37r/2 < arg(s) < ll7r/6 since, in the limit |s| 00, the 
convergence of ( |l2|) requires Re{s^) < 0. The waving 
curves are the branch cuts of the integrand in (O) . 



ImS 



x<0 



x>0 



ImS 



Saddle Point 



Branch Cut 



Saddle Point 



ReS 



SDC 




Saddle Point 
CI 



^ Saddle Point 

FIG. 9. Contour of Laplace integral satisfying B.C. Left 
figure is for a: < 0. Right figure is for a; > 0. 

The SDC of saddle point P/x approximates the con- 
tribution from the contour C3. We can find, however, 
that this contribuiton is obtained by integral represen- 
tations of some special functions more correctly than by 
the SDC approxiamtion. So we use SDC approximation 
only for the contributions from the contours Ci and C2. 
The dominan t de pendence on x of these contributions are 
obtained by ( |A^ ) to give us (pi]). 



To carry out the contour integral of C3, we transform 
s to t by sa; = e^'^^t. While we choose the argument 
of s as — TT < arg(s) < tt, the argument of t is to be 
< arg(i) < 2tt. Then (|l|) becomes 



03 (a;) 



-in—illauj^ /12k 



-P/2+iaui^/2K a/2-1 



exp 



dti-t)-^'/^ (1 + ^) 

\ —IbJX J 



a 
2k 



3x3 



'2^ 



X 



(A5) 



where Ct is the contour on f-plane corresponding to C3, 



and iuj = e^^^^uj and 



-ILOX 



J3-!v/2 



ljjx. We can simplify 



the exponent of ( |A5| ) to be exp(— <) in our parameter 
region (|9|) by the same argument in [Q: by 1/|kx3| <c 1 
and u)x <^ 1, the correction by exp[(a/2K)(- • •)] ~ 1 + 
{a/2K,){- ■ •) is negligible. With dividing the contour Ct 
into Ca and Cb as shown in the figure ^ and using 
the integral representation of Whittaker function, we can 
express (A5) as the sum of two Whittaker functions. 




Im t 



CA 



Re t 



Re t 



CB 

FIG. 10. Contour of Laplace integral for 03 (a;). 

Further we can simplify it by one confluent hypergeo- 
metric function to be 



e-"^/^ 2 sinh 

r 



^(1 



—i— 1 

K V 2 



X 



F 



1 



(A6) 



where ^ — ioLUS^ j^K^ /3 = 1 + iuj/k, T is the gamma 
function and F is the confluent hypergeometric function. 
Then we obtain the dominant dependence on x of (/)3, 
which bocomes (p2|). 
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